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SUMMARY 

Two-dimensional calculations were made of spark-ignited premixed-charge 
combustion and direct-injection stratified-charge combustion in gasoline fueled 
piston engines. Results are obtained using kinetic-controlled combustion sub- 
model governed by a four-step global chemical reaction or a hybrid laminar 
kinetics/mixing-controlled combustion submodel that accounts for laminar kinet- 
ics and turbulent mixing effects. In this work, the numerical solutions are 
obtained by using KIVA-II computer code. The KIVA-II code uses a kinetic- 
controlled combustion submodel governed by a four-step global chemical reac- 
tion; (i.e., it assumes that the mixing time is smaller than the chemistry). 

Our efforts involved the implementation of a hybrid laminar/mixing-controlled 
combustion submodel of Abraham et al. into KIVA-II. In this model, chemical 
species approach their thermodynamics equilibrium with a rate that is a combi- 
nation of the turbulent-mixing time and the chemical-kinetics time. The com- 
bination is formed in such a way that the longer of the two times has more 
influence on the conversion rate and the energy release. An additional element 
of the model is that the laminar-flame kinetics strongly influence the early 
flame development following ignition. 

Future efforts involve the implementation of this hybrid combustion model 
to study the combustion processes in High Speed Civil Transport engines and 
small gas turbine engines and to better predict N0 X emission and unburned 
hydrocarbon . 


INTRODUCTION 

The overall objective of this research is to implement a state-of-the-art 
laminar and turbulent combustion model into a three-dimensional compressible 
Navier-Stokes code. The KIVA-II code has been chosen for use in the present 
work since it is the most developed of the available multi-dimensional com- 
bustion computer programs for application of the in-cylinder combustion dynam- 
ics of internal combustion engines. There are features of KIVA-II that make it 
well suited for other applications as well, for example gas turbine combustors. 
By starting from an existing code KIVA-II and implementing an existing hybrid 
laminar kinetics/mixing-controlled combustion model of Abraham et al . (ref. 1) 
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into KIVA-II and assessing the accuracy of this hybrid combustion model by 
comparisons between predictions obtained from kinetic-control led and hybrid 
laminar/mixing-controlled combustion models, results can be obtained at much 
earlier dates. Furthermore, this hybrid combustion model has already been 
applied successfully to predict spark-ignited spray combustion in a rotary 
engine (refs. 1 to 3) and premixed-charge combustion in reciprocating engines 
(ref. 4). In this work, predictions are obtained for both spark-ignited 
premixed-charge and Di rect- Inject ion St rat i f ied-Charge combustion in gasoline 
fueled piston engines. The results of this study will be utilized in future 
works which will involve the implementation of this hybrid combustion model and 
its extension to include multi-step chemical reactions and chemical species to 
study the combustion processes in gas turbine combustors. This is part of an 
effort to develop/use advanced computer models to analyze and design combustor 
components and subcomponents, understand the physics and determine how to opti- 
mize the design to improve the performance of High Speed Civil Transport engine 
and small gas turbine engine combustors. 

Recent advances in high-speed and large-memory computers, numerical algo- 
rithms, grid generation, and mathematical modelling of turbulent reactive flows 
have markedly expanded the horizon of turbulent reactive flows that can be 
studied by numerical methods. However, full numerical simulation (or direct 
simulation of turbulence) and large-eddy simulation of a simple turbulent flame 
model still require a great deal of computer time and memory. Thus, it is 
necessary to introduce models into the codes for processes that occur on time 
and length scales that are too short to be resolved on the grid and time step 
used by the code. The use of models to describe unresolved physical processes 
introduces some empiricism into the calculations and the validity of the models 
can be assessed by comparison with experimental data. 

Combustion processes and modelling has typically been treated by means of 
(1) complete reaction-rate models, (2) eddy-breakup and mixing-controlled 
models, (3) equilibrium chemistry, and (4) finite-rate chemistry (or kinetic- 
controlled) models. In complete reaction models (ref. 5), it is assumed that 
instantaneous reaction occurs at any region where both fuel and air are pres- 
ent; the extent of the reaction is determined by the stoichiometric limit, 
i.e., the amount of fuel present determines the extent of the reaction under 
fuel-lean condi t ions, whereas the amount of oxygen determines the extent of 
the reaction under fuel-rich conditions. Complete reaction models are mixing- 
controlled and assume that the mixing time is much longer than the charac- 
teristic chemical kinetics time. This assumption may not be valid in the 
recirculation zones, and near solid walls such as combustor walls. However, if 
the flow is decelerated and the pressure and temperature increase so that reac- 
tions can proceed faster than mixing processes, a flame may be formed which is 
controlled by mixing and diffusion. 

Equilibrium chemistry models (ref. 6) reduce the number of conservation 
equations for the chemical species to be solved. It is assumed that the chemi- 
cal kinetics time is zero and reactions are mixing-controlled so that as soon 
as fuel and oxidizer are present at the same region, reactions occur as long as 
the mixture is within its flammability limits. 

In eddy-breakup models, the reaction rate is determined by the entrainment 
of fuel and oxidizer into the flame front. This entrainment depends on the 
local levels of turbulence, and the mean concentrations of fuel, oxidizer and 
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products, or the root-mean-square of the concentration fluctuations (refs . 7 
and 8). One of the main drawbacks of eddy-breakup models is their inability to 
predict local chemical equilibrium. Local chemical equilibrium can be achieved 
by expressing the reaction rate as the ratio between a concentration difference 
and the mixing (or chemical) time (ref. 9), The concentration difference is to 
be interpreted as the difference between the local concentration and that cor- 
responding to local chemical equilibrium at the same local temperature. 

Finite-rate chemistry models assume that reaction rate is kinetics- 
controlled. These models assume that the mixing time is much smaller than the 
chemical kinetics time. These combustion models are not adequate when turbu- 
lent mixing is slower than the chemistry; however, they can be used satisfac- 
torily to predict chemically controlled reactions (i.e., low temperature 
reactions, quenching, etc.). These models can be classified in two groups, 
global or quasiglobal and detailed chemistry models (refs. 10 and 11). Global 
models reduce the number of kinetic equations to be solved. However, these 
models may be able to predict qualitatively the average value of heat release 
and unable to predict local chemical equilibrium unless the reaction rates are 
expanded in Taylor series expansion around their local equilibrium values. 
Detailed chemical kinetics mechanisms are accurately known only for simple 
hydrocarbons, e.g., methane (ref. 12), propane (refs. 10, 11, and 13), and for 
hydrogen-oxygen combustion (ref. 14). For simple hydrocarbons, these mecha- 
nisms involve about several hundred chemical reactions and numerous chemical 
species, and have only been used to predict simple flows (e.g., complete 
stirred reactors) (ref 13), and one-dimensional flame propagation phenomena. 
These studies (e.g., well-stirred reactors computations) do not account for 
turbulence-chemistry effects and do not provide the flow field information. 
Nguyen and Ying used detailed propane chemical kinetics mechanisms to perform 
benchmark and sensitivity studies of spray combustion of two-dimensional gas 
turbine combustion (ref. 11). 

Depending on the nature of the time scales involved in the chemical reac- 
tions occurring in different regions inside a combustor, the relative roles of 
chemical kinetics and mixing in turbulent combustion can be interpreted by the 
above combustion models for limiting cases of very fast or very slow reactions, 
where the chemical source terms and their relation to turbulent mixing can be 
modelled. In the present study, a hybrid laminar kinetics/turbulent mixing 
model of Abraham et al . (ref. 1) is implemented into KIVA-II, and this combus- 
tion model is tested on both premixed-charge combustion and direct-injection 
s t rat i f i ed-charge combustion in gasoline fueled piston engines. This hybrid 
combustion model selects the reaction rate within the combustion flow field 
region based on a comparison between the local chemical kinetics and mixing 
times. In gas turbine combustors and internal combustion engines, combustion 
would occur in both the kinetic- and mixing-controlled regimes in different 
combustor regions; therefore hybrid combustion models would be required to 
accurately predict the combustion flow field. In the current study, results 
obtained using hybrid combustion model are compared with predictions obtained 
using kinet ic-control led combustion model. In this work, the hybrid kinetics/ 
mixing combustion model is limited to quasiglobal reaction mechanisms. How- 
ever, detailed chemistry models would be included in future works. 
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PROBLEM FORMULATION 


Engine descriptions . - The specifications and test conditions for the 
piston-cylinder configuration are listed in table I. For these calculations, 
the numerical grid consists of 20 grid points in the axial direction and 
22 grid points in the radial direction. Both the direct-injection stratified 
charge and premixed-charge engine were tested at an overall equivalence ratio 
of about 1.0. 

The calculations were made using the KIVA-II computer program, which 
solves three-dimensional equations of transient turbulent chemically reactive 
fluid dynamics with spray using an Implicit-Continuous Eulerian (ICE) technique 
for flow solver. The equations and the numerical method of solution are dis- 
cussed in detail by Amsden et al. (ref. 15). Turbulence is modelled using the 
two-equation k-e model equations. The standard turbulence model constants 
Ci, C2, and C3 (ref. 15) were used and their values are given in table II. A 
Bessel function velocity profile (ref. 15) was used to simulate the swirl 
velocity profile of the incoming air. Standard boundary conditions were used 
in this two-dimensional axisymmetric engine model and the details are described 
in reference 15. For the direct-injection strati fied-charge gasoline engine 
cases, liquid fuel (isoctane) was injected starting at crankshaft angle 5.2° 
BTDC for an injection duration of 12.6°. 

Details of the hybrid kinetics/mixing combustion model and modification 
of KIVA-II to implement this combustion model are reviewed in the following 
section. 


Combustion models . - The hybrid kinetics-mixing controlled combustion 
model of Abraham et al. (ref. 1) was incorporated into KIVA-II. This combus- 
tion model simulates the growth of the initial flame kernel into a fully devel- 
oped turbulent flame and its subsequent propagation. The reaction rates are 
proportional to the difference between the local concentrations and their cor- 
responding equilibrium values and inversely proportional to a characteristic 
time required for the achievement of equilibrium. With this hybrid kinetics- 
mixing combustion model, the reaction rates are controlled by chemical kinetics 
(or turbulent mixing) if the characteristic mixing (or chemical kinetics) time 
is much shorter than the chemical kinetics (or mixing) time. 

In the species conservation equation the time rate of change of species 
mass fraction due to conversion from one species to another is given by 


dY. 
1 

dt 



(l) 


where Y. is the mass fraction of species i and Y? is the local and 

instantaneous thermodynamic equilibrium value of the corresponding mass frac- 
tion. The present hybrid kinetics-mixing combustion model considered isooctane 
oxidation including the formation of carbon monoxide and carbon dioxide. Seven 
chemical species were considered in the hybrid combustion model: CsHi8> O2 - 

^2 . H2O, CO2 , CO, and H2 . Nitric oxide was also accounted for in the present 
calculations. N0 X was assumed to be present only in trace amounts and was 
treated spearately from the energy-releasing species. The N0 X model used was 
based on the extended Zeldovich mechanism (ref. 16). This model assumes that 
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the NO formation rate depends on the temperature, the local equilibrium compo- 
sition of the C-O-H system, and the local NO concentration. Fuel-bound nitro- 
gen N0 X is not considered in this model. Future works will attempt to extend 
the hybrid combustion model to incorporate more detailed reaction mechanisms 
including N0 X reactions and additional chemical species. 

The characteristic time x c appearing in equation (1) for the achievement 
of equilibrium is assumed to be the sum of the local laminar kinetics time and 
the local turbulent mixing time. This ensures that the rates of conversion of 
reactants to products of the reactions considered are controlled by the longer 
of the local turbulent mixing or local laminar kinetics time 

x = x n + fx . (2 ) 

c 2. t 


The expressions for x^ and x^ are 


and 


if 


if 


where 



x 


t 



k 

e 



( 3 ) 


( 4 ) 


( 5 ) 

( 6 ) 


( 7 ) 
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The delay coefficient in equation (2) is equal to unity except during ignition, 
as will be discussed below. The constants of the laminar kinetics time equa- 
tion (3), the function (1 + C|<j> - 1 . 15 1 ) and the constants B’ , C, and n 
were obtained by calibrating computed and measured laminar flame speed over a 
range of equivalence ratios (4> ) , pressures (P), and temperatures (T). The 
kinetic constants used in this work are for isooctane (ref. 4) and they are 
given in Table II. Tq and Pq are reference temperature (298 K) and pressure 
(l atm), respectively, 4> denotes equivalence ratio. One laminar kinetics time 
was used for all the seven chemical species considered. 

The turbulent mixing time in equations (4) and (6) is the same as that of 
Mangnussen and Hjertager (ref. 17), modified for the formulation of equa- 
tion (1). The constants C and C are given in table II. The subscript 

m 2 m 3 

P denotes a sum over all of the products (CO, CO 2 , H 2 , and H 2 O), and the sub- 
script PS denotes the product sum at the time of the spark (ref. 4). Equa- 
tions (4) and (6) used to describe the turbulent mixing time are generally 
applicable to diffusion as well as to premixed flames, and the local turbulent 
mixing time is controlled by the smaller of the fuel, oxygen, or products local 
concentrations. The subscripts F and O 2 denote fuel and oxygen, respec- 
tively. The superscript * denotes local thermodynamic equilibrium value. 

The delay coefficient f appearing in equation (2) has been given the 
following form 


f = 1 - exp 



and 


x d " 


c a 

m i 


(8) 


(9) 


(t - t§) is the time after spark, and the delay coefficient f sets the ini- 
tial value of the conversion time x c equal to tjj, (laminar kinetics time) and 
delay complete consideration of turbulent mixing until the time required for 
the laminar flame to traverse a distance equal to C times the length scale 
fi. of the turbulent eddies (3), ra l 


2 . 


c 3/ 4 

V 


e 


where = 0.09 (3). The laminar flame speed was given as (18) 


(10) 


S L ’ B m - B 2<» - V 2 ^ 1 - 2A r > 


(ID 


where a = 2.18 - 0.8(<J> - 1), 3 = -0.16 + 0.22(4> - 1), r is the residual mass 
fraction, and the constants B ra , <J> m , and B 2 are 0.263, 1.13, and 0.847 for 
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isooctane (ref. 4). The delay time appearing in equation (8) accounts for 

spark ignition effect on laminar flame (ref. 19). 

For this work, results obtained using chemical kinetics-controlled combus- 
tion model are calculated using the four kinetics reaction mechanisms (includ- 
ing the extended N0 X Zeldovich mechanism and six equilibrium reactions. These 
reactions and species are described in reference 20. Also, a fast algebraic 
solver, CHMQGM (ref. 15), was used to evaluate the local equilibrium species 
concentrations for the same equilibrium reactons (ref. 20) used in both calcu- 
lations employed either hybrid combustion model or kinetics-controlled combus- 
tion model . 


RESULTS AND DISCUSSION 

Comparisons of combustion rate in a direct-injection stratif ied-charge 
gasoline fueled piston engine computed using kinetics-controlled combustion 
model (case 1) and hybrid combustion model (case 2) were made at various crank- 
shaft angles (TDC is at 0). The comparisons were also made for a case involv- 
ing hybrid combustion modelling of a premixed charge gasoline engine with same 
engine dimension (case 3). In all the cases studied, the simulations were 
started at the beginning of compression stroke and ended at the crankshaft 
angle of 65° ATDC. 

Figure 1 shows the computed mean gas velocity, turbulent kinetic energy, 
and fuel spray at crankshaft angle about 29° BTDC for cases 1 to 3. The veloc- 
ity vectors are similar for all three cases and are predominantly in the direc- 
tion of the piston motion. Figure 1(b) shows the turbulent kinetic energy in 
the cylinder at the same crankshaft angle. Under each contour plot are the 
high (H) and low (L) values of the contours and the intervals of the lines that 
are between the H- and the L-lines. It may be seen that the turbulent kine- 
tic energy is high for all the engine cases studied with the piston cup with a 
bowl. The added azimuthal swirl profile also enhances the turbulent kinetic 
energy. This would increase the vaporization, fuel-air mixing and burning 
rates. Figures 1(a) and (c) indicate that the regions of high turbulence lev- 
els increase the fuel vaporization rate. At crankshaft angle about 11° BTDC, 

4° BTDC, and about TDC, the increase in density increases the tubulent kinetic 
energy as TDC is approached (see figs. 2(b), 3(b), and 4(b)). This would 
increase the burning rate and flame propagation as TDC is approached; and this 
effect is more pronounced for case 2 (hybrid combustion model). Comparisons 
between cases 1 and 2 of the turbulent kinetic energy profiles, the fuel vapor 
distributions and gas temperature profiles for each of these crankshaft angles 
indicate that the hybrid combustion model predicts faster burning rate in 
regions of high turbulent kinetic energy and temperature, and therefore less 
amounts of unburned fuel vapor are present in these regions. It would be 
expected that the turbulent mixing time (tj) would control the species conver- 
sion rates and energy release in these regions using the hybrid combustion 
model. As also seen in these figures, the regions marked by low tubulent kine- 
tic energy (e.g., cylinder head and wall) are responsible for the poor mixing, 
and results in incomplete combustion as observed by high levels of unburned 
fuel vapor. These figures also show that the temperature profiles and unburned 
fuel vapor profiles near the cylinder head and wall are similar for case 1 
(kinetics-controlled) and case 2 (hybrid combustion model). This is caused by 
the kinetics character isic time (tj,) controlling the species conversion rate 
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as the flame approaching the low temperature regions (cylinder head and wall), 
and similar predictions of the burning rate would be obtained using either the 
kinetics-controlled combustion model or the hybrid combustion model. 

Some simplified observations can be made for the premixed-charge combus- 
tion engine flowfield. Premixed combustion results in rapid and high energy 
release (temperature rise) as seen in figure 2 (case 3) and more uniform burn- 
ing. The burning rate (as observed by the isotherms) can be correlated with 
the spatial distribution of the turbulence in the engine. High turbulence 
intensity results in faster burning rate as observed by the hot flame kernel in 
the regions of high turbulent kinetic energy. As the premixed turbulent flame 
approaches the cylinder wall, piston face, and piston bowl, the decrease in 
turbulence intensity slows down the flame propagation resulting in slow and 
incomplete combustion in these regions (see figs. 2(b), (d), and (e), 3(b), 

(c), and (e), and 4(b), (d), and (e) for case 3). 

At crankshaft angle about 15° ATDC, most of the liquid fuel has vaporized 
(cases 1 and 2) and most of the fuel (cases 1 to 3) has been depleted. 


CONCLUSIONS 

Kinetics-controlled combustion model and hybrid kinetics-mixing controlled 
combustion model have been used to make comparisons of mixing and burning rate. 
The comparisons have been made for a direct-injection strati f ied-charge combus- 
tion and a premixed-charge combustion gasoline fueled piston engines. The 
hybrid combustion model incorporates relaxation effects due to chemical kine- 
tics and turbulent mixing times. The results obtained for the direct-injection 
strati f ied-charge combustion engine indicate that burning rate predicted by the 
hybrid combustion model is faster in regions of high turbulent kinetic energy, 
resulting in high temperature flame kernel and more depletion of the fuel. As 
the flame approaches the cylinder head and wall, the laminar kinetics time 
would eventually become controlling, thereby similar burning rate is obtained 
for the direct-injection strat i f ied-charge engine using the kinetics- 
controlled or hybrid combustion model. 

The hybrid combustion model with the same model constants has also been 
tested for a premixed-charge combustion engine and yielded reasonable results. 

Future works will extend this hybrid combustion model to include more 
chemical species and the implementation of this combustion model into KIVA-II 
to study the combustion processes in High Speed Civil Transport and small gas 
turbine engines and to more accurately predict N0 X emission and unburned 
hydrocarbon. 
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TABLE I. - ENGINE PARAMETERS AND TEST CONDITIONS 



DISC engine 

Premi xed-charge 
combustion engine 

Engine stroke, cm 

9.55 

9.55 

Compression ratio 

6.54 

6 . 54 

Volume of piston bowl, cm^ 

58.8 

58.8 

Squish clearance, cm 

0.1819 

0.1819 

Connecting rod length, cm 

16.269 

16.269 

rpm 

1600 

1600 

Initial engine air temperature, K 

400 

360 

Initial engine pressure, atm 

1 

1 

Fuel 

Isooctane 

Isooctane 

Starting crankangle of injection 

52° BTDC 


Injection duration, deg 

12.672 


Injection angle, deg 

60 


Starting crankangle of ignition 

27° BTDC 


Ignition duration, deg 

9.6 



TABLE II. - MODEL CONSTANTS 


Turbulence: C^ = 

1.44* 

*"2 = 1.92; C 3 = 1 

,*-io 

Combustion model : 

c 

= 2.75; C = 0.055; 

C = 0.092; B' = 4.59x10 

m l 

m„ 



2 

3 


i = 

1.0; n = -0.75; E /R 

3 

= 15098 K; C = 0.08 
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CASE 2 


CASE 3 


CASE 1 



VELOCITY ACROSS J = 1 PLANE, -29.50°, 
U MAX = 7-84828E+02, V MAX = 2.94169E+03 
W MA v = 9. 76695E+02 



V VvO\ '■ \ \ ' 
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VELOCITY ACROSS J = 1 PLANE, -28.98°, 
U max = 7 .83083E+02, V MAX = 2.94321E+03 

W max = 9 . 74576E+02 

(a) GAS VELOCITY. 


I I , 



' ■ i J / / - " ‘ 



VELOCITY ACROSS J = 1 PLANE, -26.31°, 

= 8.21862E+02, V MAX = 2.98933E+03, 

W MAX = 9. 15953E+02 



■ TKE ACROSS J = 1 PLANE, L = 4.61711E+04, 
H = 2.98336E+05, 29.05°, DO = 3.15206E+04 



.TKE ACROSS J = 1 PLANE, L = 4.60790E+04 
H = 2.97494E+05, 28.98°, DQ = 3.14269E+04 



.TKE ACROSS J = 1 PLANE, L * 4.57243E+04, 
H - 2.93886E+05, DQ = 3.10203E+04 


(b) TURBULENT KINETIC ENERGY. 



‘1 



305 PARTICLES IN THE SYSTEM 


316 PARTICLES IN THE SYSTEM 


(c) FUEL SPRAY. 

FIGURE 1. - RESULTS AT CRANKSHAFT ANGLE ABOUT 29° BTDC FOR CASE 1 (KINETICS-CONTROLLED COMBUSTION MODEL; DISC ENGINE), CASE 2 
(HYBRID COMBUSTION MODEL; DISC ENGINE) AND CASE 3 (HYBRID COMBUSTION MODEL; PREMIXED CHARGE PISTON ENGINE). 
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VELOCITY ACROSS J = 1 PLANE, -11.37°, 
U MAX = ® ■ 84606E+02 , V MAX * 3.38367E+03, 

W MAV = 1 . 24565E+03 
max 



VELOCITY ACROSS J = 1 PLANE, -13.45°, 
U MAX = 8.01494E+02, V MAX = 3.23956E+03, 

W MAV = 1 . 33204E+03 

MAX 


(a) GAS VELOCITY. 



VELOCITY ACROSS J = 1 PLANE, -12.60°, 
U MAX = 9. 11390E+02, V^ = 2.72781E+03, 

W MA¥ = 8.86937E+02 

MAX 



.TKE ACROSS J = 1 PLANE, L = 4.80499E+04, 
H = 2.39197E+05, -11.37°, DO = 2.38934E+04 


.TKE ACROSS J = 1 PLANE, .TKE ACROSS J = 1 PLANE, L - 6.56371EKJ4, 

L = 5 . 38007E+04, H - 2.73481E+05, H = 3.05680E+05, -12.60°, DO = 3.00053E+04 

-13.45°, DQ = 2.746006+04 

(b) TURBULENT KINETIC ENERGY. 




122 PARTICLES IN THE SYSTEM 133 PARTICLES IN THE SYSTEM 

(C) FUEL SPRAY. 

FIGURE 2. - RESULTS AT CRANKSHAFT ANGLE ABOUT 11° BTDC FOR CASES 1, 2, AND 3. 
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c 8 h 18 ACR0SS J = 1 PLANE, L * 1 . 17529E-02, 
H = 1.05776E-01, -11.37°, DO = 1.17529E-02 



.TEMP ACROSS J = 1 PLANE, L = 6.57971E+02, 
H = 2.96074E+03, -11.37°, DO = 2.30277E+02 


CASE 2 



c g h 18 ACROSS J = 1 PLANE, L = 9.20729E-03, 
H = 8.28655E-02, -13.45°, DQ = 9.20729E-03 


(d) FUEL VAPOR. 



.TEMP ACROSS J = 1 PLANE, L = 7.127535+02, 
H = 2.64682E+03, -13.45°, DQ = 1.93407E+02 

(e) GAS TEMPERATURE. 

FIGURE 2. - CONCLUDED. 


CASE 3 



C g h 1g ACROSS J = 1 PLANE, L = 6.59493E+03, 
H = 5.93543E+02, 150, -12.60°. 

DQ - 6.59492E-03 



.TEMP ACROSS J = 1 PLANE, L = 7.37012E+02, 
H = 2.90261E+03, -12.60°, DQ = 2.16560E+02 


OF POOR QUALITY 





CASE 1 


CASE 2 


CASE 3 





VELOCITY ACROSS J = 1 PLANE, -4.60°, 

= 1 .021 12E+03, = 3 . 24035E+03 , 

H hax = 1*98914E+03 





VELOCITY ACROSS J - 1 PLANE. -5.80°, 
U MAX = 9.57021E+O2, V MAX = 3.08009E+03, 
W MA¥ = 1.06227E+03 

MAX 

(a) GAS VELOCITY. 



VELOCITY ACROSS J = 1 PLANE, -3.92°, 
U MAX = l *’W*1E+03, V^x = 6.35114E-K)3, 

W kax = 3.06003E+03 



.TKE ACROSS J = 1 PLANE, L = 5.13525E+04, 
H = 3.3897SE+05, -4.60°, DQ * 3.59529E+04 


.TKE ACROSS J = 1 PLANE, L = 6.40585E+04, 
H = 3.91696E+05, -5.80°, DO = 4.09547E+04 

(b) TURBULENT KINETIC ENERGY. 


.TKE ACROSS J = 1 PLANE, L = 2.39673E+05, 
H = 1.92218E+06, -3.92°, DQ = 2.10313E+05 




82 PARTICLES IN THE SYSTEM 90 PARTICLES IN THE SYSTEM 

(C) FUEL SPRAY. 

FIGURE 3. - RESULTS AT CRANKSHAFT ANGLE ABOUT 4° BTDC FOR CASES 1, 2 AND 3. 
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CASE 1 


CASE 2 


CASE 3 



u 

m! 

L 


Cghig ACROSS 
H = 8.32176E 



J = 1 PLANE, L = 9.24640E-03, 
02, -4.60°, DQ = 9.24640E-03 


c ghig ACROSS J - 1 PLANE, L - 6.18474E-03, Cgh lg ACR0SS J = 1 PLANE, L 

H = 5.56626E-02, -5.80°, DQ = 6.18474E-03 H = 1.08816E-03, -3.92°, DQ 

(d) FUEL VAPOR. 



.TEMP ACROSS J = 1 PLANE, L = 7.91687E+02, .TEMP ACROSS J = 1 PLANE, L * 

H = 2.76721E+02, -5.80°, DQ * 1.97553E+02 H = 3.32848E+03, -3.92°, DQ = 

<e) GAS TEMPERATURE. 



C0 2 ACROSS J = 1 PLANE, L = 9.68919E-03, co 2 ACROSS J = 1 PLANE, L = 7 

H = 8.72027E-02, -5.80°, DQ - 9.88919E-03 H = 6.93867E-02, -3.92°, DQ - 

(f) co 2 . 

FIGURE 3. - CONCLUDED. 


= 1.20906E-04. 
= 1.20906E-04 


2.63391E+03, 

6.94568E+03 


. 70964E-03, 
7.70964-03 
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CASE 1 
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CASE 3 



VELOCITY ACROSS J = 1 PLANE, 3.33° 


Vx = 1-091UE+O3. V MAX 
= 1 . 28796E+03 


= 2.85898E+03, 
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VELOCITY ACROSS J = 1 PLANE, 0.98°, 

U MAX = 1.03939E+03, V^ = 3.20517E+03, 

W MAY = 8 . 29506E+02 

MAX 

(a) GAS VELOCITY. 
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VELOCITY ACROSS J = 1 PLANE, 4.56°, 

U MAX = 1 .24574E+03, V^ = 4.06404E+03, 
V MX = 2.32013E+O3 



.TKE ACROSS J = 1 PLANE, L = 6.59326E+04, 
H = 4.69689E+05, 3.33°, DQ = 5.04696E+04 



.TKE ACROSS J = 1 PLANE, L = 5.93149E+04, .TKE ACROSS J = 1 PLANE, L = 1.60359E+05, 

H * 4.04443E+05, 0.98°, DQ = 4.31410E+04 H = 1.323S0E+06, 4.56°, DQ = 1.45405E+05 

(b) TURBULENT KINETIC ENERGY. 



65 PARTICLES IN THE SYSTEM 

FIGURE 4. 


53 PARTICLES IN THE SYSTEM 
(C) FUEL SPRAY. 

RESULTS AT CRANKSHAFT ANGLE ABOUT 3° ATDC FOR CASES 1, 2 AND 3. 
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CASE 1 


CASE 2 


CASE 3 




c 8 h 18 ACROSS J = 1 PLANE, L = 5.12489E-03, 
H = 4.62591E-02, 3.33°, DQ = 5.14178E-03 



.TEMP ACROSS J = 1 PLANE, L = 7.95260E+02, 
H = 3.00501E+03, 3.33°, DQ = 2.20976E+02 


c 8 h lg ACROSS J = 1 PLANE, L = 5.23541E-03, 
H = 4.71 187E-02, 0.98°, DQ = 5,23541-03 

(d> FUEL VAPOR. 



.TEMP ACROSS J = 1 PLANE, L = 8.33670E+02, 
H = 2.80571E+O3, 0.98°, DQ = 1.97204E+02 

(e) GAS TEMPERATURE. 


c 8 h 18 ACR0SS J = 1 PLANE, L = 1.76187E-07, 
H - 1.58568E-06, 4,56°, DQ = 1.76187E-07 


'if. 


.TEMP ACROSS J = 1 PLANE, L = 2.49589E+03, 
H = 3.20936E+03, 4.56°, DQ - 7.13467E+01 


/ 


co 2 ACROSS J = 1 PLANE, L = 1.04958E-02, co 2 ACROSS J = 1 PLANE, L = 7. 4347 IE -03, 

H = 9.44623E-02, 0.98°, DQ = 1.04958E-02 H - 6.69124E-02, 4.56°, DQ = 7.43471E-03 

(f) co 2 . 

FIGURE 4. - CONCLUDED. 
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VELOCITY ACROSS J = 1 PLANE, 15.66°, 
U hax = 1 .2M358E+03, - 2.64965E+03, 



VELOCITY ACROSS J = 1 PLANE, 12.72°, 
U MAX = 1 . 38928E+03 , V^ = 3.00436E+03, 

= 8.82616E+02 



VELOCITY ACROSS J = 1 PLANE, 19.93°, 
U MAX = 1 .21416E+03, U NAX = 3 . 13883E+03, 

W max = 1 -G61 I »0E-K>3 



.TKE ACROSS J = 1 PLANE, L = 6.21923E+OM, 
H = *I.M51*l2E+05, 15.66°, DO * M.78687E+OM 


.TKE ACROSS J = 1 PLANE, L = 5.27897E+04, 
H = 3. 7251 ME +05, 12.72°, DO = 3.99655E+04 


.TKE ACROSS J = 1 PLANE, L = 1.03728E+05, 
H = 8. 16490E+05, 19.93°, DG = 8.90952E+04 


(b) TURBULENT KINETIC ENERGY. 




48 PARTICLES IN THE SYSTEM 


2 PARTILES IN THE SYSTEM 


(c) FUEL SPRAY. 

FIGURE 5. - RESULTS AT CRANKSHAFT ANGLE ABOUT 15° ATDC FOR CASES 1, 2 AND 3. 
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CASE 1 


CASE 2 


CASE 3 



Cgh 1g ACROSS J = 1 PLANE, L = 5.74846E-14, Cgh^g ACROSS J = 1 PLANE, L = 3.77290E-03, Cgh 1g ACROSS J = 1 PLANE, L = 2.67012E-03, 

H = 5. 17576E-13, 19.93°, DQ = 5.75114E-14 H = 3.39584E-02, 15.66°, DQ = 3.77319E-03 H = 2.40307E-02, 12.72°, DO = 2.67008E-03 

(d> GAS TEMPERATURE. 



.TEMP ACROSS J = 1 PLANE, L = 7.61224E+02, .TEMP ACROSS J = 1 PLANE, L * 7.88991E+02, .TEMP ACROSS J = 1 PLANE, L = 2.15037E+O3, 

H = 2.82021E+03, 15.66°, DQ = 2.05899E+02 H = 2.67605E+03, 12.72° DQ - 1.88706E+02 H = 2.89086E+03, 19.93°, DQ = 7.40485E+01 

(e) GAS TEMPERATURE. 



CO 2 ACROSS J = 1 PLANE, L = 1.11053E-02, co 2 ACROSS J = 1 PLANE, L = 1.03770E-02, 

H = 9.99480E-02, 12.72°, DQ = 1.11053E-02 H = 9.33934-02, 19.93°, DQ = 1.03770E-02 

Cf) co 2 . 

FIGURE 5. - CONCLUDED. 
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